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A two-dimensional nonrelativistic fermion system coupled to both electromagnetic gauge 
fields and Chern-Simons gauge fields is analysed. Polarization tensors relevant in the quan- 
tum Hall efi^ect and anyon superconductivity are obtained as simple closed integrals and are 
evaluated numerically for all momenta and frequencies. The correction to the energy density 
is evaluated in the random phase approximation (RPA), by summing an infinite series of 
ring diagrams. It is found that the correction has significant dependence on the particle 
number density. 

In the context of anyon superconductivity, the energy density relative to the mean field 
value is minimized at a hole concentration per lattice plaquette (0.05~0.06)-(pc''/^)'^ where 
Pc and a are the momentum cutoff and lattice constant, respectively. At the minimum 
the correction is about —5 % ~ —25 %, depending on the ratio 2mwc/p^ where lj^ is the 
frequency cutofi^. 

In the Jain-Fradkin-Lopez picture of the fractional quantum Hall effect the RPA correc- 
tion to the energy density is very large. It diverges logarithmically as the cutoff is removed, 
implying that corrections beyond RPA become important at large momentum and frequency. 



1. Introduction 

There is a growing interest in Chern-Simons theory. Gauge fields with a pure Chern- 
Simons term, which we call Chern-Simons gauge fields, effectively alter the statistics of 
matter fields. [1-5] In Jain's picture of the fractional quantum Hall effect, [6] an electron 
and magnetic flux form a bound state, or a composite fermion. The fractional quantum 
Hall effect for electrons is understood as integral quantum Hall effect for these composite 
fermions. Further, Fradkin and Lopez have shown that these composite fermions are re- 
garded as fermions interacting with Chern-Simons gauge fields. [7] Halperin, Lee, and Read 
have advocated this picture further to investigate physics near half filling in the quantum 
Hall effect. [8,9] 

In anyon superconductivity, physics of essentially the same system of fermions coupled 
to Chern-Simons and Maxwell gauge fields is explored, but with a different Chern-Simons 
coefficient. [10-26] It is known that this system shows a rather unique behavior at finite 
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temperature, particularly around T ^ lOOK. [25,26] So far there has been no evidence for 
its possible connection to high Tc superconductors. In fact, the theory of anyon supercon- 
ductivity is still in its infancy to the extent that it provides only few predictions which can 
be tested experimentally. 

In the relativistic case the Chern-Simons gauge fields induce boson-fermion transmuta- 
tion. [27] If gauge fields have both Maxwell and Chern-Simons terms, the Chern-Simons term 
gives the gauge fields a topological mass. [28] Further it has been recently shown that the 
presence of Dirac fermions leads, under certain conditions, to the spontaneous breakdown of 
the Lorentz invariance by dynamical generation of a magnetic field (B ^ 0).[29] There quan- 
tum fluctuations play a crucial role in lowering the energy density of the true ground state 
with B 0. Lopez, Rojo, and Ftadkin have shown that an effective relativistic Dirac theory 
with Maxwell-Chern-Simons gauge interactions naturally arises in the quantum Heisenberg 
antiferromagnet on a square lattice. [30] 

In this paper we examine a nonrelativistic fermion system interacting with both Maxwell 
and Chern-Simons gauge fields. Employing the technique developed in ref. [29], we first 
cast one-loop polarization tensors for gauge fields in simple integral forms suited for both 
analytical and numerical evaluation. Then we derive the exact formula for the energy density 
in terms of gauge field propagators. Applying our one-loop formula for the polarization 
tensors, which corresponds to the random phase approximation (RPA), or to summing 
up an infinite series of ring diagrams in perturbation theory, we examine contributions of 
quantum fluctuations to the energy density. 

The detailed numerical evaluation is given in the absence of dynamical electromagnetic 
interactions. We shall find that RPA corrections generate non-trivial dependence of the 
energy density on the particle (or hole) number density in the context of anyon supercon- 
ductivity. The energy is minimized at a number density precisely where high Tc material 
is superconducting. In the case of the fractional quantum Hall effect, RPA gives a large 
correction to the energy density, in accord with the equivalence between the electron picture 
and the composite fermion picture. 

2. The model 

The model we consider consists of a nonrelativistic fermion field ip, electromagnetic field 
A^, and Chern-Simons gauge field a^. Its Lagrangian is given by 

£ = i {eE] - xB') - ^ e'^'^'a^d^ap + i^^Do^ - ^{Dki>)^ {Dki>) + epAo 

Ej = Foj , B = -Fi2 = diA^ - d2A^ (^.1) 
Dn = df,+ i{afj, + eA^) . 

The last term represents a contribution from a background neutralizing charge. Note that 
a'^ = — ttfe and A'' = —Ak. The Euler equations are 

edkEk = e{f-p) 
e'^^xdiB =e/ 

-Is^'^Pf - (2-2) 
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where the currents are given by 

f = v-V 

i (2.3) 

j'' = -—{i^^Dki>-iDk^)H} . 
2m 

Suppose a uniform external magnetic field -Bext is applied and the system remains trans- 
lation invariant. Eq. (2.2) implies that ip feels a total magnetic field 

6tot = 6^°) + eSW = ^ + eSext • (2.4) 

K 

The Landau level densities associated with the external magnetic field and the total magnetic 
field are given by 

„ext _ |e-Bext| 

Pl - 



1,^'; (2.5) 

„tot _ rtoti 
Pl 



2tt 

Something special happens when the total filling factor is an integer: 

^=T^=" • (2-6) 
Pl 

The change in the statistics phase induced by Chern-Simons gauge fields is AOs = 1/2k. 
Combining (2.4), (2.5), and (2.6), one finds 

In the Jain-Fradkin-Lopez picture, the system exhibits the fractional quantum Hall effect 
when AGs is a multiple of 27r. Suppose that eSext > 0. Then the condition is satisfied if 

^ = -2. 

_ ^ (2.8) 
2pn±l 

where ± corresponds to e(6tot)- In the case eSext < the signs of k and btot are reversed. 
The main sequence in fractional quantum Hall effect is given by p = 1, or z/cxt = n/ (2n± 1). 
In this paper we focus on the special case v = 1. Generalization to the case = n is 
straightforward, but is left for a future investigation. Note that for v — 1 



e{btot)}p=-^ . (2.9) 



3. Fermion propagator 

The perturbation theory is defined around 6tot- We rewrite a^°^ + and A^^ — > 

A^'^' + A^, i.e. and A^ now represent fiuctuation parts. Then the Lagrangian becomes 

^ = -f ^ext + -Cr'l^, A] + rr^'^lV', V-l + jCintia, A, i,, ^t] . (3.1) 
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£gauge _|_ ^matter (jefines the zeroth order "free" part. The gauge field part is given by 

£gauge ^ 1 (^^2 -^B')-^ e^'^Pa^d^ap + C^,. [a, A] (3.2) 

where £g.f . [a, A] is a gauge fixing term to be specified in the next section. The matter field 
part is 

The interaction part is given by 

-Cint = - xBe^tB - (V'V - P) (ao + e^o) 

The BextB term in Cint, being a total derivative, does not contribute in perturbation theory 
and may be dropped. 

The zeroth order fermion propagator has been well described in the literature. Usually 
it is given in the form of an infinite sum over the Landau level index. Our task here is to 
recast it in an integral form for later convenience. 

Let us quote the result from Ref. [26]. In the Landau gauge 

1 _ (3.5) 

p - \^tot\ ■ 

I is the magnetic length. When v lowest Landau levels are completely filled, the fermion 
propagator is given by 

G{x,y) = -i{T[^ix)i^\y)]) 



Go{x) = - J 0(.xo) -e{-xo) J2 \ e""""" 

^ n—iy n=0 

1 f°° 

^ J dze-"'''/^ Vn[z-z{x2)]v„[z + z{x2)] 



(3.6) 



where 

1\ 1 



en=(n+-)^ (n = 0,l,2,...) (3.7) 
is the n-th Landau energy level, and 

_ (-1)" c.V2_^ 

^"W 2"/27ri/4(n!)i/2'' dx^" (38^ 
z{x2) = e(6tot) |y • 

The phase factor in the propagator G{x,y) is not translation invariant, but does not con- 
tribute to any physical quantities as we see below. 
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We now utilize the Feynman-Hibbs summation formula[31] 

~ Qit/2 fir 1 1 

V e-*"* Vn{x)vniy) = , . . exp <^ — - (x^ + y^) cost - 2xy \ (3.9) 

in order to perform the infinite sum and obtain an explicit integral expression for the prop- 
agator. We also quote another formula 

y'±^!:hM=, r dt'-j==L exp / [(^2 + cost- 2x2/1 1 (3.10) 
^n-A-^e Jo \/2msint '^i2sina^ ' J ' 

which follows from (3.9). 

First we evaluate the propagator in the v = Q case. Employing (3.9) and taking the 
Fourier transform 

Go(p) = Go{uj,p) = Jdx Go{x)e'^^-°-P'^^ , 

one finds that 

Goip)\u=0 = -imp e-i(pl+pl)l't.n{t/2)+i^mlH _ (311^ 

Jo cos(t/2) 

For i/ > 1 we observe that 

Go{x)\. - Go{x)\,=o 

'-'-1 ■ <.oo 

Using the explicit expression for u„, one finds that for i/ = 1, for instance, 
Go(p)U=i =47rie-p"'' ^(w - eo) 

-imZ^ '^^ ^-»p tan(f /2)+ic^m;^t 

io cos(i/2) 



4. Gauge field propagators 

Free gauge field propagators are obtained from Cf^^^^ in (3.2). For the sake of unified 
treatment of both electromagnetic and Chern-Simons gauge fields, let us consider gauge 
fields described by 

Afo = l {eE^ - xB^) - \ e^^PA^d^A, - ^ (£^0^° + xd^A^^f ■ (4-1) 

The last term represents gauge-fixing. Electromagnetic fields are described with k = 0, 
whereas Chern-Simons fields in the e, x ^ limit: 

£r''k^]=M)[A;K = 0]+A^[a;e,x-0] . (4.2) 

As is obvious, the propagator is singular in the e, x ^ limit, which is merely a gauge 
artifact. Physical quantities such as the energy density are well defined in the e, x ^ 
limit. We take this limit when convenient and safe to do so. 



(3.12) 
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In passing, we mention that for pure Chern-Simons gauge fields the alternative gauge 
fixing term 

>Cg.f. = ^ {dka'r (4.3) 

is also very useful and convenient. [26] The radiation gauge is obtained in the a — > limit. 
We adopt the gauge fixing in (4.1) in this paper. 
The Lagraiigian (4.1) has the form 

/{e/x)V Kd2 -KdA 

K^"" = ~Kd2 ~V Kdo , ^^-^^ 

\ ndi —kOq —V J 
D = e9o'-xV' . 

For this Lagrangian, the propagator, the inverse of K^^" , is straightforwardly found to be 






Xd2 


-Xdi 


-Xd2 





edo 




edo 






^^°^- = ^PT^j|^2 7 + « e^o 1^.(4.5) 

Equal time commutation relations are given by 

[A^(t,x),n'^(t,y)]=iV^(x-y) 

[A^{t,x),A,{t,y)] = [n''(t,x),n-(i,y)] = 



X 



W = edoA, - ^e^^Ak 

In evaluating the energy density we need the full gauge field propagators. The interaction, 
given by (3.4), mixes the electromagnetic and Chern-Simons gauge fields. The propagators 
have a matrix structure: 

When the lowest Landau level is completely filled so that the perturbative ground state is 
of Hartree-Fock type, one can develop a Feynman diagram method. Since both gauge fields 
couple to fermions in the combination of + eA^, one can write 

_(D^^-,^ \_ (e- e\ (4-8) 

Here T^'^ represents the sum of one-particle irreducible diagrams common to both and 
eAa. 
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Gauge invariancc, translational invariancc, and rotational invariance imply that T^'^ is 
expressed in terms of throe independent invariant funetions life's: 

+ (1 - - -p'<5^")(n2 - Ho). ^ ■ ' 

All life's are functions of and p^. 



5. Evaluation of Ilk's 

We evaluate the kernel r^''(w, g), or equivalently the invariant functions life's, to the 
leading order in the case v = \. The interaction is given by (3.4). One-loop contributions 
are depicted in Fig. 1. 

To calculate life's we take a frame q = {q, 0). Then 

Diagram (a) is easily evaluated to give 

m 

As is well known, other diagrams yield 
T^'^°°{q)=il-0^Go{p)Go{p-q) 

r^'^'' (9) = / {DtG,{p) ■ D^Goip - 5) + D+Goip) ' DtGo{p - q) ^5-2) 

+ D+D^Goip) ■ Goip -q) + Go{p) ■ D^D+Go{p - q)] 

D^Goip) = {iP2 ± ^(^tot)^^) Goip) ■ 

The phase factor in G{x,y) (3.6), does not contribute to F'"^. 

The next step is to insert (3.13) into (5.2) and perform the momentum integrals. All 
integrals involve heavy, but similar manipulations. We shall describe the computation of fig 
in detail. 

There are four terms in the integral 

I- Jo cosp i 

X (47rie-(P^-«)''' 5{po - uj - eo) - iml^ H ^-i{p-,n'i'tMt'/2)+iipo-^)mi^t''[ _ 

I- ^0 cos it' J 

(5.3) 
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The product of the two J-function terms is easily evaluated to be 

-ie-«^'V2^(^) . (5.4) 

Two cross terms of the (5-function and t-integral pieces are, after po integration, 

-i{p-q)^l^ ta.n{t/2)+i{eQ-u)ml^t 



J (27r)3 Jo 



dt _,Vi?_;r^2;2 J 



COS 



e 



.m f°° dt 1 , ,2 ig2|2tank\ ,^ ^, 

/ r: TT:t- IT exp ^(eo - ^^)mlH - — — — j- (5.5) 

^TT Jo COS 2* 1 + * tan 2^ V l + ztan2r/ 

+ {-uj,-q)} 

- r dte-'''''^'-^-''y' cosiivmlH) . 
The product of the two f-integrals in (5.3) vanishes: 

Jo cos U Jo cos k' J 



ni 
= I- 



(27r) 



j dtdt' J dp ■■■6{t + 1') ■ 



(5.6) 







The manipulation in the last equality is justified since the integrand is regular &tt = t' — 

no(g) is a regular function of u so that the 5{(jj) term in (5.4) must be cancelled by a 
part of (5.5). To see this we rewrite (5.5) as 



m 



POO 

Jo 



_ m r 

~ 27rio 



/o 5t 1^ (w + ie)mP (w — ie)ml 

Integration by parts yields 



2nP \lj 



hi) 



1 f 1 

ie LO — ic 

277 Jo V 2 J I (t^ + if)mP [ijj - ie)mP 

Now the second 5{u)) piece is 



2 /'oc 2 roo 

^S{uj) (iie-^*-^*-«'''(i-'="")/2 = %5(a;) / (-ids) 
2 Jo 2 Jo 



e 



s-g^i-'(l-e-')/2 



(5.8) 



^(e-^'V2_i)^(^) . 
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Hence (5.5) becomes 

i e-«'''/2 5{uj) r e-^*-«'''(i-«"")/2 sm{ujmlH) . (5.9) 

The first term exactly cancels (5.4). Combining (5.4), (5.6), and (5.9), one finds that 



1-K 



Ho = — ^b(|?^/^a;mZ^ 



(5.10) 

Fo{x,y) = / rft sinyte-"-^(i-'^"") . 

y Jo 

The evaluation of IIi and 112 proceeds similarly. The cancellation of J-function pieces 
takes place as in the case of Ho. The result is 

""^'^^"'^'^"''''"""'^ 

Q (5-11) 
Fi{x,y) = —x-Fo{x,y) 



F2{x,y) = -l + !^x+(^-^xy^Fo{x,y) . 



In refs. [11,13,21,26,32] these life's are obtained in the form of infinite sums over the 
Landau level index n. Applying the second Feynman-Hibbs formula (3.10) to the result in 
ref. [26], the results (5.10) and (5.11) are reproduced with some labor. 

Notice that all life's are related to one function Fq. For small x [33] 



oo u 



Fo(x,y) = e-^Y^^ • (5.12) 

fc! t/2 _ (fc + 1)2 + ie ^ ^ 



In particular 



Fo{0,y) 



1 — y^ — ie 



Foix,0) = — [ 

X Jo 



dw — 



(5.13) 



w 

One can derive an alternative integral representation for Fj which is more suited for 
both analytical and numerical evaluation. First we Wick-rotate the y variable: Fj {x, z) = 
Fj{x,iz). Then 

Fo{x,z) = ■ h{x,z) 

°° ^k+i ^ (5.14) 

h{x,z) = 



h{x,z) satisfies 



x-^^ + z"^ ) h{x,z) = xe^ 



h{Q,z) = Q (5.15) 
dxh{Q, z) = . 
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In terms of s = In a;, h{x, z) = J*^ du p{u) smz{s — u) where p{s) = xe^. Or 



hix, z) = - dw sin (z In — ) • e"" 
z Jo \ wJ 

X 

= — / dt sin(2lnt) -e^* . 

z Jo 



(5.16) 



This h{x,z) satisfies the boundary condition in (5.15). 
Hence one finds that 



-t)x 



1 

Fo{x,z) = — / dt 8m{z\nt) ■ e~^^~ 
z Jo 

1 /•! 

Fi{x,z) = — / (it{l + x(t- 1)} sin(zlni) -6^(1-*)^ (5.17) 
^ Jo 

1 

l + F2{x,z) = — / dt{l-2x + 2,xt + x^{t-lf} sm.{zhit)-e~ 

z Jo 



-(l-t)x 



The behavior of these functions is depicted in Fig. 2. The asymptotic behavior is given by 

^ 1 



, ^ X 1 Sx^z^ 

1 + F2 1 . 

x"^ + x"^ + z"^ (x^ + z^)^ 

These approximate expressions are vahd to the accuracy of 2 % for x or z > 10. 
6. Energy density 

Quantum fluctuations shift the energy density of the ground state from the mean flcld 
value. This shift can be related to the fuU gauge field propagators. We first derive an exact 
formula for the energy density [34], and estimate it by utilizing the result in the previous 
section. The mean field energy density for 1/ = 1 is given by 



^mean = ^^^ext + ' (6-1) 
^ TTv 

Here Bext is related to p and k by (2.9). The interactions (3.4) give corrections to fmean- 
The corresponding interaction Hamiltonian is 

= Jd^[{a'' + eA'')^(i;Wk^|J-{Dk0i^)+{ao + eAo){^P^^|J-p)] 2) 

if(2) = [ dx-^(a'= + eA'=)2^V • 
J ^th 

Now we consider a Hamiltonian given by 



H{g) = Ho+gHW+g^H(''^ 



(6.3) 
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where Hq is the free part of the Hamihonian. With an auxihary parameter g, H{g) connects 
the mean field Hamiltonian -ff(O) and the full Hamiltonian H(l). Its ground state satisfies 

H{g)\^{g)) = E{g)\^{g)) 

E{g) = {^{g)\H{g)\^{g)) . 

The desired change in energy density with inclusion of the interactions is E{1) — E{0). Since 
(*(5)|*(5)) = 1, 



Eil)-EiO)=l\g±EU 



j\g{n9)\^^m9)) (6.5) 



= f dg{^{g)\H(^^ + 2gH^^'>\^{g)) . 
Jo 



The currents associated with the fluctuation parts of the gauge fields are defined by 



a. = r 



e,X=0 



(6.6) 



■c=0 

Ki^" has been defined in (4.4). In the theory descibed by the Hamiltonian II{g), 
? = 5(^t^-p) 

Zm m 
Straightforward substitution of these currents yields 

j dx j>(x){a^(a;) + eA^{x)] ^ gH^^^ + 2g^H^^^ . 

Thus 



(6.8) 



E{1) - E{0) = /' ^ / dx ( *(ff)| ?(a^ + eA^) ^9) ) • (6.9) 
Jo 9 J 

If the gauge fields A^ have (4.1) as a free Lagrangian, their equal- time commutation 
relations are given by (4.6). Elementary algebra shows that 

KPi^{T[A^{x)A,{y)] ) = {T[KP'^A^{x)A4y)] ) 

+ ^^''{l + ^^0(^ - l) ]eSixo - yo) [A^{x), A^y)] (6.10) 
= (T[K'"^A^{x)A,{y)] ) + iS'\S\x - y) . 
The difference between the full and bare propagators, therefore, satisfies 

Kp/[{T[A^{x)AM) - {'^[M^)My)])o} = {'^[K''^M^)My)]) ■ (e.n) 
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Applying (6.6) and (6.11) to (6.9), one obtains that 
E{1) - E{0) 

= J^'^Jd^ lim (^K^'-''^^^{{T[A,{x)A.{y)] ) - {T [A^{x)A.{y)] )„} 
+ i^cS'-''{(T [a^(ar)a.(j/)] ) - ( T KWa.(y)] )o} 

Fourier-transforming this expression yields, for the density, 

£il)-£iO) 



(6.12) 



= 'l'jj ^{^'D^,em{D[AA] - D^^) + trDo-^5(i?H - 



(6.13) 



. dg f d^p 



9 



Tr r- 



Dq-i - f 



Here the two-by-two matrices D, Dq, and f arc given by (4.7) and (4.8). The trace tr 
is taken over only Lorentz indices, whereas Tr is taken over both Lorentz and gauge field 
species. We stress that the formula (6.13) is exact, no approximation being involved. 

We have evaluated F^^ to the leading order in Section 5. With the Hamiltonian (6.3), all 
Ilfc's there must be multiplied by g"^. In other words, F = 0{g^). The g-integral in (6.13), 
then, can be easily performed: 



(6.14) 



F^^^ is precisely the kernel evaluated in Section 5, which represents the sum of one- loop 
diagrams. The final expression in the above equation has simple diagramatic interpretation. 
It represents the sum of a series of ring diagrams connected by either Chern-Simons or 
electromagnetic fields. (See Fig. 3.) The 'In' takes care of combinatoric factors. 
In the frame = {u, q, 0) 



g^IIo LuqUa iqlli 
V'^ip) = I Loqila lu^Ho iujUi 
-iqlli — zwlli 112 



(6.15) 



With a general propagator Dq given by (4.5), 

/ -Xg2(no + -Hi) ea;g(no + -ni) iq{eUi + KUo)\ 



-TDo 



es — 



-Xwg(no + ^Hi) eto^ilLo + ^ILi) iLo{eIli + kHq) 
^ ixq{Ui + ^Ti2) -iew(ni -I- ^na) eUi + KUi ^ 



(6.16) 
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where s = euj"^ — x'f' ■ The electromagnetic part is obtained by taking the k ^ Umit in the 
above formula. For the Chern-Simons part we take the limit e, x ^ with the ratio /? = e/x 
fixed. We find 



iqbi —iujb2 d 

ai = Ho Hi , 02 = — IIo III (f\i7\ 

e s KU s KU y^-'-') 

Xe^TT ^TT f. ^^TT 1 TT 

Oi = III 112 , 02 = — Hi 112 

e S KU S KU 

1 e^- 1 

c=— ni--no , rf = — n2--ni , 

S K S K 

where u = ui^ — j3q^. The /3-dependence in the expression for the Chern-Simons part should 
disappear in physical quantities. 

Indeed, for the energy density (6.14), we have 



Af = £{1) - £{Q) 

d^p 



The shift in the energy density in pure Chern-Simons or pure Maxwell theory is found as a 
special case: 



Af|pu.e cs = -| / ^ In ^ {(Hi - Kf - Uon^} 

|p..e BM = 4 / ^ In { (i + ^ no) (i + ^ n2) - ^ n? 

In terms of Fj{x,z) = Fj{x,iz) introduced in (5.10) and (5.11), 

-=^r-r-n{i.^.2.^.o-^A 



(6.19) 



(6.20) 



where s = —{2x/l'^)x — {e/w?l'^)z'^ . Here we have introduced ultraviolet cutoffs, h.^ and A^, 
supposing that the model (2.1) is an effective theory valid at low energies. They are related 
to the momentum and frequency cutoffs by A^^ = p^Z^/2 and A^ = ml'^cOc- 

Expression (6.20) contains many parameters. Let us take anyon superconductivity as a 
typical example, p ^ 10^^ cm^^ so that 6tot ^ 100 T and / ^ lOA. The lattice spacing a 
is about 5A. For high Tg superconductors one expects Pc ~ and Uc ~ leV. The bare 
electron mass is m = 2.6 x 10^°cm~^. With these values A^; 2, and A^ ~ 2.6. Note that 
Az/Aa; = 2mu)c/Pc s-nd that the effective mass m may be different from the bare electron 
mass. The coupling constant e'^/An is ~ a/d where a = 1/137 is the fine structure constant 
and d 5A is the interplanar spacing. This gives e^/47r ^ 1.5 x lO^cm"^, e'^l/An 0.015 
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and e'^ml'^/ATT ^ 40. However, in the Fradkin-Lopcz picture of the fractional quantum Hall 
effect, there are no cutoffs. We shall come back to this point in Section 8. 
In pure Chern-Simons theory 

—2 

A5=^i?(A,,A,;c) 
m 

i?(A^,A^;c) = - / dx dzln^{x,z;c) (6.21) 
71" Jo Jo 

^{x,z;c) = {l-cFif-c^FQF2 , c-- 



\2 „2 „ _ e(^tot) 

27rK 



Here we have made use of the relation v = 1 and p = (27ri^)~^. R{Ax,Az;c) represents 
the RPA correction relative to the mean field value. It has an important dependence on 
the number density p through A^; and A^. A detailed numerical study is presented in the 
following sections. 

The behavior of R{Ax, A^] c) at large A^; and A^ is analytically estimated with the aid 
of (5.18). One finds that for A^,, A^ > Aq > 10 



R{Ax, A^; c) =i?(Ao, Aq; c) + tan ^ ^) 



c2 



|-ln^ — — / du — tan ^ u\ . 
[2 Ao TT Ji u J 

Notice that it grows logarithmically as A, or as the density p gets smaller. 



(6.22) 



7. Anyon superconductivity 

A neutral anyon gas with pure Chern-Simons interactions has an application to anyon 

superconductivity. [10-26] Previously, the energy density in the neutral anyon gas has been 
evaluated in the Hartree-Fock approximation by Hanna, Laughlin and Fetter. [12,35] The 
correction to the energy density is — ^ (— ^) of the mean field value in the ly = 1 [i/ — 2) 
case, independent of the particle number density p. 

It is known that high Tc ciipratc superconductors become superconducting only when a 
hole concentration Xc, namely the number of holes per lattice plaquette, is in the range 0.05 
to 0.25. Within the context of anyon superconductivity it has been an unanswered question 
why superconductivity is achieved only in a limited range of Xc- 

We shall show that in RPA the energy density is minimized exactly at a number density 
in the range mentioned above. At a lower or higher number density the energy is increased, 
and the system is expected to lose superconductivity. 

Since c = 1 for z/ = 1 and Bext = 0, the total energy density of a neutral anyon gas is, 
from (6.21), 

f(p)*°* = ^{l + i?(A„A,;l)} 

i?(Aa,, A^; 1) = - / dx dz \n^{x,z;l) . 
TT Jo Jo 

The density dependence comes in through A^; = p^/Anp and A^ = ma;c/27rp. R is evaluated 
numerically. 
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First wc depict, in Fig. 4, the behaviour of the argument of the logarithm, z\ 1), in 
the integrand. It vanishes at the origin, behaving as 2x + -^x^ + z^. It is a smooth function, 
approaching 1 as a; or gets large. 

As explained in the previous section, the momentum and frequency cutoffs in anyon 
superconductivity are approximately Pc ^ and Wc ~ 1 eV, where a is the lattice spacing. 
This gives A^/Aj; ^ 1, but there remains an ambiguity. We have evaluated i? as a function 
of Aj;, or equivalcntly of {pca)^'^Xc = p/Pc = l/(47rA2;), with Az/A^ = r fixed. 

Fig. 5 (a) shows the behavior of R for a wide range of Aa;. When Ax,A^ > 10 (small 
density) it grows logarithmically (~ ilnA^;) as expected from (6.22). 

In Fig. 5 (b), R is plotted in the range < Xc/iPco)'^ < 0.5 where the minimum takes 
place. We recognize that the minimum is located around xd {pcd/'hY = 0.05 ^ 0.06 for 
a wide range of r (0.2 < r < 5). \i Pca/h ~ 2, x™™ ~ 0.2. Although the value of is 
ambiguious, it is safe to conclude that at a lower number density the energy density sharply 
increases and therefore that the superconductivity is lost at low densities. The correction at 
the minimum ranges from —5% to —25%, depending on r. It is rather surprising, but also 
is encouraging, that the minimum occurs approximately at a number density where cuprate 
material is superconducting. 

Certainly, more elaboration and detailed study of the anyon superconductivity model is 
necessary in order to see if it can describe high superconductors or yet-to-be-discovered 
new material. We leave it for future investigation. 



8. Fractional quantum Hall effect 

In the Jain-Pradkin-Lopez picture of the fractional quantum Hall effect, electrons in 
an external magnetic field with a filling factor i/gxt are described by composite fermions 
interacting with, in addition to the external magnetic field, Chern-Simons gauge fields, 
as was explained in Section 2. [6-9] Our previous result is valid for the total filling factor 
u = n = 1, which applies, from (2.8), to a sequence fext = l/(2p ± !)• 

In the electron picture the lowest Landau level has an energy |eBext|/2m so that the 
mean field energy density is 

felectron IQ -\\ 

I'ext rn 

In the composite fermion picture the lowest Landau level has an energy |6tot|/2m so that 
the mean field energy for i/ = 1 is 

^composite 2^ 

Clearly f^-^^"^'*'' < E^^^°^. There appears a large discrepancy, by a factor of 

Since the transformation from the original electron system to the Chern-Simons-fermion 
system is exact, the discrepancy observed above is merely an artifact of the mean field 
approximation. We are going to show that quantum fluctuations give a big correction. It is 
expected that if one includes all corrections, the energy densities in the two systems should 
be exactly the same. 

Without loss of generality we take eSext > as before. From (2.8) one finds that for 

zy = 1 

.ext = ^ Where ± = e(6,oO 
c = T'^p . 



Page 16 



D. Wesolowski, Y. Hosotani and S. Chakravarty 



As shown by Halperin, Lcc and Read [8] and by Simon and Halpcrin [9] , RPA gives 
a substantial correction to the excitation energy spectrum. In the absence of the Coulomb 
interaction, the spectrum is determined by zeros of detD^g cx (Hi — k)^ — 110112. In terms 

(1 - cFi)2 - c^FoFa = . (8.4) 
At a zero momentum x = \q^l^ = 0, 

Fo(0,2/) = Fi(0,2/) = l+F2(0,y) 

^ 1 (8.5) 
1 — y'^ — ie 

Substitution of (8.5) into (8.4) yields 

y = toZ^o; = |c- 1| = 2p± 1 = — , (8.6) 

t'ext 

or 

a;(g = 0)=^PA = — ^ . (8.7) 

This is exactly the Landau gap in the original electron picture. The RPA correction yields 
a factor z/^j compared with the mean field value in the Chern-Simons picture. 
The energy density in RPA is given by 



I M,, .A, (8-8) 
R{h.x-,^z'-,c) = — I dx I dz \a.^{x,z;c) . 

Jo Jo 

We have evaluated i? as a function of A^; for various values of r = ^z/^x and c. There 
correspond two values of c to each fext- For instance, u^^t = 5 (^) corresponds to c = — 2 
and +4 (—4 and +6). The function lR.{x, z; c) in (8.8) is depicted in Fig. 6. 

For a large cutoff A^:, R ~ ^c-^lnA^; as follows from (6.22). Its global behaviour is 
depicted in Fig. 7. For a given i/ext, the RPA correction exceeds the saturation value u~J^ — 1 
at around Aa, = 10 100. This implies that RPA is not accurate for large momenta and 
frequencies. 

In the light of the equivalence between the original electron picture and the composite 
fermion picture the cutoffs must be removed. In one respect our result shows that higher 
order corrections beyond RPA become crucial in the computation of the energy density. It 
asserts that quantum fluctuations give a large correction of order v~J^ . 



9. Summary 

We have examined the nonrelativistic Maxwell- Chern-Simons gauge theory relevant in 
the quantum Hall effect and anyon superconductivity. We have derived compact integral 
representations for the polarization tensors r^'^{p) or Ilk{p)- Quantum fluctuations give 
important contributions to the energy density of the system, which is most conveniently 
expressed in terms of the gauge field propagators. The energy density was evaluated in 
RPA. 
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One of the important consequences in RPA is that there resuhs a non-trivial dependence 
of the energy density on the particle number density. In particular, in anyon supercon- 
ductivity, the ratio of the energy density in RPA to that in the mean field approximation 
is minimized at a hole concentration Xc = (0.05 ~ 0. 06) {p ca / h)"^ . In light of the appli- 
cation to high Tc superconductors this is extremely intriguing and encouraging. Typical 
high Tc cuprates become superconducting for 0.05 < Xc < 0.25. Neither the mean field 
approximation nor the Hartree-Fock approximation have been successful in explaining why 
superconductivity is achieved only in a limited range of Xc- 

In the application to the fractional quantum Hall effect, the RPA correction to the energy 
density in the Pradkin-Lopez picture diverges as the cutoffs are removed. The divergence 
comes from large momenta and frequencies where RPA is expected to break down. We may 
conclude that higher order quantum fluctuations give substantial corrections to the energy 
density, presumably restoring the equivalence between the original electron picture and the 
composite fermion picture. 

In this paper we have concentrated on the case of a unit filling u = 1 with respect to 
the total magnetic field 6tot = p/k + eSext- The generalization to arbitrary integer filling 
v = n is important both for the anyon superconductivity and for the fractional quantum 
Hall effect. It is left for future investigation. 
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Figure captions 

Fig. 1 One- loop corrections F'^" defined in (4.8) and (5.2). Tliree kinds of vertices are given 
by the last three terms in (3.4). 

Fig. 2 The behaviour Fj{x,z) in (5.17). (a) Fi{x,z), (b) F2{x,z) and (c) 1 + F2{x,z). 
All arc smooth functions of x and z = iy. They start out with unity at the origin 
and approach zero asymtotically for large x or z (5.18). 

Fig. 3 The energy density in RPA. See Eq. (6.14). A doubly dashed line represents the sum 
of Chern-Simons and electromagnetic field propagators. 

Fig. 4 Plot of the argument of the natural logarithm in (7.1) as a function of x and z = iy. 
It vanishes at the origin and approaches unity asymtotically for large x ot z. (See 
text for more details on the asymtotic behaviour.) This is related to the RPA energy 
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density of a neutral anyon gas. Note that conventional order of x and z has been 
switched for visual clarity. 

5 Magnitude of the RPA correction to the mean field energy density R{\xt^z] 1) for 
a neutral anyon gas normalized to the mean field value for various choices of the 
lattice cutoff A^^ and A^; (7.1). The solid line denotes cutoff ratio r = A^/A^; = -g, 
the dashed line denotes r = 1 and the dotted line denotes r = 5. (a) Plotted as 
a function of A^; for fixed values of r on a logarithmic scale. Linearity for large 
values of A^; is due to the logarithmic growth of R with cutoff (6.22). Departure 
from linearity at very large values of A^; is understood to be a numerical artifact. 
Note that the minima for r = 1 and r = 5 are not resolved. They are resolved 
in the following figure. Numerical estimates of R are accurate to about 5%. (b) 
Alternative presentation of Fig. 5(a) around the minima. R has been plotted as a 
function of Xc{pcCi/K)~'^ = l/(47rAa:) (scaled hole concentration), where Pc and a are 
the momentum cutoff and lattice spacing, respectively. The minima are adequately 
resolved and their locations are more or less independent of the cutoff ratio r. For 
large values of concentration, R approaches zero. We note intriguing similarity with 
the cuprate superconductors. Numerical estimates of R are accurate to about 10%. 

6 Plot of the argument of the natural logarithm in (8.8) for i/gxt = \ functions of 
X and z = iy. Corresponding values of the constant c are —2 and +4. Both start 
with 3? = 9 at the origin. See text for more details on the asymtotic behaviour. Note 
that conventional order of x and z has been switched for visual clarity, (a) Plot for 
c = +4. (b) Plot for c = -2. 

7 Magnitude of the RPA correction to the mean field energy density R{Kx,h.z;c) for 
a quantum Hall system normalized to the mean field value for various choices of 
the lattice cutoff A^, and A2 (8.8). R is plotted as a function of for fixed values 
of r = Aji/Aj; on a logarithmic scale. Linearity for large values of A^; is due to 
the logarithmic growth of R with cutoff (6.22). R is always positive in all cases. 
Departure from linearity at very large values of A^^ is understood to be a numerical 
artifact. Numerical estimates of R are accurate to about 5%. (a) Behaviour for a 
c = 4 (t'cxt = \) system. The solid line denotes cutoff ratio r = i, the dashed line 
denotes r = 1 and the dotted line denotes r — b. (b) Behaviour for various choices 
of external filling fraction with r = 1. The solid line denotes c = — 2 (i/ext = \)j the 
dashed line denotes c = +6 (i^ext = \)^ the dotted line denotes c = — 4 (z/ext = \) 
and the dash-dotted line denotes c = +4 (z/gxt = |)- 
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